Hybrid Luttinger liquid/Fermi liquid behaviour in an atomic wire on a surface 
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Recent advances in single atom manipulation have made it possible to create "wires" in the form 
of atomic scale linear structures on a semiconductor surface. If such structures are to be used in 
future electronic devices, it will be necessary to know to what extent they behave as one-dimensional 
conductors, in which the presence of electron-electron interactions in the wire lead to properties quite 
different from those found in higher-dimensional systems. In particular, it would be useful to know 
if these structures retain any of the Luttinger liquid properties that are predicted for a pure one- 
dimensional metal. However, experimental studies of such structures have so far yielded unclear and 
contradictory results. We investigate this problem theoretically by creating a highly simplified model 
of a wire on a surface — a chain of atoms that includes electron-electron interactions to represent 
the wire, coupled to a similar chain that excludes electron-electron interactions to represent, albeit 
crudely, the surface. We present results for the eigenstates and spectral functions for such systems 
that suggest that many Luttinger liquid indicators are present. However the system as a whole 
retains some residual Fermi liquid characteristics, and therefore merges properties of both one and 
higher dimensions. 

PACS numbers: 71.10.Pm, 73.22.Lp, 74.20.Mn, 79.60.Jv, 81.07. Vb 



I. INTRODUCTION 

In recent years it has become possible to create atomic 
scale wires consisting of a single chain of atoms on a 
semiconductor surface. These wires have obvious poten- 
tial applications in future nanoscale circuitry, but most 
of the research effort so far has concentrated on the elec- 
tronic structure of these systems. 

Moreover, it has been shown that certain of these 
systems are not subject to the Peierls transition, even 
at low temperatures. Hence they are expected to be- 
have as one-dimensional metals, with properties substan- 
tially different from the Fermi- liquid properties of higher- 
dimensional metals. Specifically, the discontinuity of the 
occupation number at the Fermi energy is lost, spin and 
charge excitations become separated and the correlation 
functions of the system develop power-law behaviour. 
Experimental evidence for such Luttinger liquid-like be- 
haviour has, however, proved elusive and ephemeral. In 
particular, the case of Au wires on vicinal Si(l 1 1) has 
attracted much recent discussion |l], g, [|, [|, ||, Q|, as 
have the potential Luttinger liquid properties of quasi- 
one-dimensional charge-transfer salts j?], ||, ||] . 

One of the important questions that remains unan- 
swered is to what extent the coupling of the wire to the 
surface affects the one-dimensional properties that would 
be expected of an isolated wire. This is a somewhat 
more complicated problem than that of two coupled one- 
dimensional conductors, which has been widely studied 
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The present work therefore begins to tackle this prob- 



lem by reducing it to its most simple form. The wire 
is represented as completely as possible, by a fully inter- 
acting one-dimensional conductor based on the Luttinger 
model. This is coupled to a one-dimensional representa- 
tion of the surface, similar to the wire but without any 
electron-electron interactions. In this way we are able to 
probe the effects of the coupling of the wire to the sur- 
face in as simple a way as possible, while retaining many 
of the properties that would be found in a more realistic 
system. 

In a previous paper (22) we described the case where 
the one-electron properties of the two chains are the 
same. In this paper we build on these results by treat- 
ing the more relevant case in which the properties of the 
chain differ. Where relevant, we have reproduced results 
from this earlier work. 

The remainder of this paper is structured as follows. 
We start with a review of Luttinger liquid theory, fol- 
lowed in section II B by a full description of our model 
and a justification of its choice. We then present results 
for the eigenstates of the system, both for two chains 
identical with respect to their non-interacting properties, 
and for two chains with differing properties. We then 
present results for calculations of the spectral functions 
for both types of systems, and calculate the Luttinger 
parameter K p . The paper concludes with a discussion of 
the implications of the results. 



II. THEORETICAL METHODS 



A. The Luttinger model 
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Our calculations are based on Haldane's solution of the 
Tomonaga-Luttinger model [ p3| , consisting of a chain of 
spinless idealized atoms with periodic boundary condi- 
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tions. The model is characterized by a linearized disper- 
sion relation with slope v-p and a complete separation of 
the populations of left- and right-moving particles (as in 
figure |lj) and the inclusion of an infinity of negative en- 
ergy electrons 23, It is this that makes the model 



exactly solvable, as the boson commutation relations are 
then exact rather than approximate. 




FIG. 1: Schematic representation of the Luttinger model. 
p = +1, — 1 indicate the populations of left- and right-moving 
fermions respectively. States up to Ef are filled, this includes 
the infinity of negative energy states. 



pq| ) and so we shall only review here those aspects of 
particular relevance to the current work. The total Lut- 
tinger model Hamiltonian can be split into three parts 



H — Hq + H 2 + Hi 



(1) 



Hq represents the non-interacting part of the Hamilto- 
nian: 



H o = ^2 v F (pk - fcp) : & pk c p k 

p.k 



(2) 



where q indicates momentum,p (= ±1) the branch index 
and : . . . : implies normal ordering. The interaction parts 
of the Hamiltonian, Hi and H4, are given respectively by 



I \ - 

H 2 = ~2^g2(q)p+(q)p-(-q), 



Hi = -^2gi(q)pp(q)p p (-q). 

P,Q 



(3) 



(4) 



The Hi term represents forward scattering between the 
left- and right-moving fcrmion branches, while the H<± 
term represents forward scattering within a momentum 
branch. The density operators p are defined in terms of 
the fermion operators as 



p P (q) = '■ £ lk+ q cp,k ■ 
k 



(5) 



The details of the Luttinger model are covered in sev- 
eral review articles (see, for example, references p4l G3, 



The complete Hamiltonian can be diagonalized via a 
Bogoliubov transformation |2j| to yield the bosonized 
form of the Hamiltonian 



J 



H = ~ J2 v ^ : Pp(q)P P (-q) ■ [vn (N+ + iV_) 2 + vj (N+ - N-f 
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u q - v F \q\) + ^u q b\b q 
1 



+ ^ (vnN 2 + vjJ 2 ) , 



(6) 
(7) 



r 



where the transformed density operators p are related to 
the originals by a phase 4> q : 



p P (q) = p P (q) cosh (j> q + p~ p (q) sinh^g 



(8) 



negative branches, and the difference between them 
(analogous to current). 

The three velocities in the Hamiltonian are related as 
follows: 



characteristic frequency 



There is 

vT(^F + 34( < z)/2 7r ) 2 — {92(q)/^) 2 \/\q\ associated with 
the transformed bosons, while the quantum numbers 
N+ — N_ represent respectively 



N 



N + + N_ and J 



the sum of the number of electrons on the positive and 



V N VJ 
V N 



v 0> 

v K p 



v e 



-24, 



(9) 



v e 



-2,p 
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and are also related to the non-interacting Fermi velocity 
by 



V N = V F + 
V.J = Vp + 



.94 + 92 

2n • 

94 - 92 

2tt ' 



(10) 



The requirement that H be a diagonalized version of H 
is ensured by the relationship between the Bogoliubov 
transformation phase 4>„ and the interaction functions 



/ 7to f + 9i{q) -32(g) 



K p (q). 



9i{q) + 92(g) 



1/2 



(11) 

(12) 



The spinless Luttinger liquid parameter K p is then ob- 
tained by taking the limit of K p (q) as q tends to zero. In 
all of the above the limit q — > is implied where q is not 
explicitly included. 

In addition there are two further parameters related to 
K p : for a spinless Luttinger liquid they are given by 



1 

7 = 2a. 



K a + — -2 

^P 



(13) 
(14) 



The Fermi liquid corresponds to K p = 1 and a = and 
so departures from these values can be used to "mea- 
sure" the extent of non- Fermi liquid behaviour, a is 
an exponent which governs the power-law dependence 
of all single-particle properties (for example the density 
of states, which varies as N(E) w \E — E-p\ a ), as well as 
other properties of the system, including the d.c. resis- 
tivity, which varies as p{T) w T 1 " 7 24 



B. The model 

We are now able to couple two Luttinger chains (la- 
belled by superscript A and B) together by allowing hop- 
ping between adjacent points on each chain with matrix 
element t±: 

^coupled -\- 



(15) 

Other schemes of inter-chain hopping are of course pos- 
sible, and we expect that they will generate similar re- 
sults. We neglect any "drag effects" of inter-chain in- 
teractions that would be generated by interaction terms 
analogous to H2 or H4 involving electrons on both chains. 
We have the choice of using either a set of basis states 
generated by the diagonalized boson operators w q or the 
non-diagonalized operators a q . We choose the latter, as 



although this has the consequence that the ground states 
for a single interacting chain no longer consist of the rele- 
vant zero-boson basis states, it is considerably more con- 
venient computationally. The ijv are fermion creation 
operators given in the bosonized form by 

1 



7? 



Jpk F x 



(16) 



z (x) is a boson field operator given by 



<t>p,c{ X ) 



1/2 



where the subscripts p refers to the branch index (±1) 
and c to the chain index. V ViC is a ladder operator whose 
form ensures the anticommutation properties of the final 
fermion operator ^(x) despite the commutation proper- 
ties of the constituent boson operators a q . 

In order for J7 P , C to produce anticommuting field op- 
erators on different chains, it is necessary to introduce a 
further phase factor into its definition, analogous to Hal- 
dane's phase factor Q(jp, N p , N^ p ) for ensuring anticom- 
mutation between the branches of a single chain. The 
total ladder operator component of equation ( |l6| ) thus 
takes the form 

U p , c = C(p, N p , AT_ P )C'( C) N c , iV_ c )|7V p , c + 1, 7V_ P , C ), 

where the subscript c = ±1 is a chain index. The anti- 
commutation properties originate in the phase factors £, 
which can be written as 

Q= P ,c = (-i)^ N - ) . (19) 

This model provides us with a substrate that is metal- 
lic, whereas in reality we want to calculate for an insu- 
lating substrate. This may be achieved within the limi- 
tations of the Luttinger model in two ways. The first of 
these, reported here, is to adjust the Fermi velocity of the 
"substrate" chain so that the electrons preferentially re- 
side on the "wire" chain. Our method of achieving this is 
described in detail in section III A . The second approach, 



to be reported elsewhere, is to move all the single-electron 
states associated with the substrate chain up in energy. 
This therefore has the effect of simulating a gap between 
the valence and (largely inaccessible) conduction bands 
of the substrate chain, while the states belonging to the 
wire chain are located within this gap. 



C. Computational details 

The choice of a specific form for the interactions 32(9) 
and gi(q) is arbitrary, as the Luttinger model remains 
completely solvable for any interaction that fulfils cer- 
tain conditions |2^]. We choose a Gaussian form for the 
interaction and set 32(9) = 9i{g) = 2itV(q) with 



V(q) = Iex V (-2q 2 /r) 



(20) 
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Periodic boundary conditions on ring of length L 



Chain A 

(with interactions) 



Interchain hopping t± between adjacent sites 



***** 



Chain B 
no interactions) 



FIG. 2: Schematic representation of the two chain model. 
Electron-electron interactions are included on just one of the 
chains. 



as this has the advantage that it maintains the same form 
in both real and momentum space. Other forms, such as 
a screened Coulomb interaction, would also be possible. 
The parameters I and r can be varied to control the 
"strength" and range of the interaction respectively. We 
complete the basic two chain model by setting the inter- 
action strength on chain B, representing the substrate, 
to zero but retaining interactions on chain A. 

The size of the computational Hilbert space is re- 
stricted by allowing only one boson in each mode. This 
does not affect the low energy properties in which we are 
interested, as we have checked that only higher energy ex- 
citations involve basis states with more than one boson 
in each mode. We also neglect the presence of the infin- 
ity of negative energy electrons, i.e. c|iV e i 0C t r ons = 0) = 0. 
However, we are still strongly limited in the size of the 
system we can handle: the Hilbert space scales exponen- 
tially with both length of chains and number of electrons. 
The largest system for which we have been able to cal- 
culate with arbitrary numbers of electrons is one with 
length L = 6 (i.e. 12 sites in total). 



D. Spectral functions 



We define the spectral function p(q, uS) as 



-Q[G^ c (k F + q,u; + p)], (21) 



with q = fep — k, and the retarded Green's function Gp c 
defined as the double Fourier transform of 



G£ c (s,t) = -i0(t){(Vw=(M)< c (o,o)>} 



(22) 



where the subscripts p and c again indicate branch and 
chain indices. Whereas for a Fermi liquid we would 
expect p(0, oj) to be a delta function at the Fermi en- 
ergy, the situation for a Luttinger liquid is quite different 
pq , p9[ |3C| . Spectral weight is repelled from the Fermi 
surface owing to the virtual particle-hole excitations gen- 
erated by the inter-branch scattering term 172, resulting 
in a broadened peak. As q is increased the Fermi liquid 
spectral function will merely broaden like q 2 reflecting 
the finite lifetime of electrons away from Ep, but for a 
Luttinger liquid there is zero spectral weight within a 



range ±voq of the Fermi energy. In addition, the neg- 
ative frequency contribution is suppressed exponentially 
with q, and for a continuum Luttinger liquid the positive 
frequency contributions have a power law dependence. 
For u> > this is of the form 6{u> — voq)(ui — unf?) 7 " 1 and 
for w < it is of the form 8(—m — vaq)(—LU — fn?) 7 , with 
7 given by equation (Jl4|) |2J] . 

However, for our more complex system it is not easy 
to obtain an analytic form and so we resort to computa- 
tional methods. We calculate the Green's function 



G R (k,k' : Lu) = (N\c k 



oj — H + en 



(N\c 



(23) 

LU + H-e N c k \N) 



using Haydock's tridiagonal Lanczos-based procedure 
pl[ |32] |. In order to ensure convergence, we put lu — > 
to + irj, where 77 is an imaginary component of the energy 
roughly equal to the level spacing of the system |n| . \N) 
is the A^-electron ground state, with energy en- We also 
use the Lanczos method to calculate the eigenstates of 
the system. 



III. RESULTS 

A. Two types of system 

Ideally, we would like to be able to calculate for a 
metallic chain coupled to an insulating chain represent- 
ing the substrate. However, the metallic nature of our 
chains is inherent in the Luttinger model and so this 
is not possible with the present system. Nonetheless, 
we can reproduce some of the most important charac- 
teristics of such a system within the constraints of the 
current model. In particular, we arc interested in sys- 
tems where the Fermi velocities of the chains differ, i.e. 
the Fermi velocity for the interacting chain is less than 



that of the non-interacting chain, 



< vi 



The re- 



sult of this is that the electrons find it energetically 
preferable to reside on the interacting chain. This is 
one of the characteristics we would expect for a metal- 
semiconductor chain system, where the electrons in the 
wire are more mobile than those in the substrate. This 
has several consequences for the band structure of the 
coupled system, as can be seen more clearly in figure 
||. The bands are not parallel and are not separated 
by 2t±, as they are for two chains with equal vf The 
bands anticross at E = 0, which is Ep for our refer- 
ence state of no participating electrons. Also, for similar 



chains, Vp — Vp = v F on mg : 
this case we have Vp(= 1.0) 



antibond 



mg 



whereas in 



v 



bonding / 



vf{= 2.0) 



V 



antibonding / 



(= 1.002) and 



(= 1.998). 



We now present our numerical results for these two 
systems — the first with Vp = Vp, the second with Vp < 
Vp. The quantities we calculate are those that may 
give us information about the Luttinger liquid proper- 
ties of the system — specifically, the boson population of 
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(a) 



Antibonding band * 
v F = 1.998/ 



/ Bonding band 
,' v F = 1.D02 




(b) 



FIG. 3: (a) Schematic band structure of an atomic wire on 
a semiconducting surface, (b) Band structure as realised in 
our calculations for the non-interacting system with L = 4, 
N = 2, t± = 0.1 Vp = 1.0, Vp = 2.0. Filled circles represent 
occupied states, unfilled circles empty states. 



the eigenstates, the excitation energies and eigenstates, 
and the spectral function. 



B. Boson contributions to eigenstates 

Using the methods detailed above, we calculate the 
ground and neutral first excited states for our two sys- 
tems at various interaction strengths. We first consider 
the boson population of these eigenstates. As we have 
chosen to work in the space of un-Bogoliubov trans- 
formed bosons, it is only the non- interacting (I = 0) 
systems that have a zero-boson population in all modes 
for the ground state. All the interacting systems (/ > 0) 
have contributions from basis states that include bosons 
on one or both chains, as these basis states are not 
individually eigenstates of the interacting system. For 
both our two chain systems, contributions to the ground 
state may come from basis states with no bosons on ei- 
ther chain, bosons on the interacting chain (chain A) or 
bosons on both chains. These contributions are plotted 



in figure |] as a proportion of the total ground state. Note 
that contributions from basis states with bosons on the 
non-interacting chain only are negligible and are not in- 
cluded. 

For the case where bot h cha ins have equal Fermi ve- 
locities, as shown in figure 4(a), the data for basis states 
with bosons only on chain A maps almost exactly onto 
that for a single isolated chain. However, there is also a 
significant contribution from states with bosons on both 
chains which increases with / — the presence of interac- 
tions on chain A induces interaction effects on the oth- 
erwise non- interacting chain B. For stronger interchain 
coupling t±, we see less of a contribution for both types, 
but nonetheless the data for bosonic states on chain A 
is still comparable to that of a single chain, indicating 
that for both t± = 0.1 and t± = 0.5 we are in a weakly 
coupled r egim e. 

Figure 4(b) compares the results for the equal vf sys- 
tem with those from the system with Vp < Vp. Once 
more, the data for bosonic basis states only on chain A 
matches the single chain data for both of our 2-chain 
systems. We again observe induced chain B boson states 
at a level similar to that for the equal vp system in the 
Vp < Vp system. 



C. Charge distribution 

We can also examine the distribution of charge be- 
tween the two chains. Taking a simple case where there 
are just two electrons available for excitation, the pos- 
sible distributions are both electrons on the interacting 
chain, one electron on each chain, or both electrons on 
the non-interacting chain. Those basis states where both 
electrons reside on the same chain may be thought of 
as "ionic" -like, and those where the electrons are shared 
between the chains as "covalent"-like. The charge dis- 
tributions for the gr ound and first neutral excited states 
are shown in figures 5(a) and 5(b). We see that for the 
ground states of the equal Vp system, the eigenstates 
evolve smoothly from the non-interacting limit (where 
the charge distribution is has equal proportions of cova- 
lent and ionic states) to a situation at large values of / 
where the ground state is dominated by chain B ionic 
states. The presence of interactions on chain A has thus 
forced all the available charge onto the non-interacting 
chain B. The case of the excited states is, however, some- 
what different. The I — excited state is a fourfold de- 
generate combination of ionic and covalent states. Once 
even small chain A interactions are switched on, this de- 
generacy is lifted and the dominant basis states are cova- 
lent, as indicated by the discontinuity at / = in figure 



5(a) . As I is increased, the role of the covalent basis 
states increases, but there is a further discontinuity in 
the charge distribution at / w 2.0, where the dominance 
of the covalent states is replaced by that of chain B ionic 
states. This is owing to a change in the nature of the ex- 
citation, from an inter-band excitation with Ak = to an 
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FIG. 4: Contribution of bosonic basis states to the ground 
states as a function of interaction strength I. 



FIG. 5: Charge distribution within the chains for the ground 
state and neutral first excited state of the two-electron, two- 
chain model. 



intra-band excitation with A/c = 2ttvf/L. For / > 2.0, 
the charge distribution for the excited state is very sim- 
ilar to that of the ground state, indicating that no net 
interchain hopping is involved in these excitations. 

For the Vp < Vp system, as shown in figure EM the 
non-interacting ground state is dominated by chain A 
ionic state, with lesser contributions from covalent states 
and a minimal contribution from chain B ionic states. 
As I is increased, these evolve smoothly through a re- 
gion (0.5 < I < 1.0) dominated by covalent states, to a 
situation, for I > 1.0, where the system is dominated, 
like the equal vf case, by chain B ionic states. The ex- 
cited states exhibit a similar pattern of discontinuities to 
the equal vf system, although the transition from inter- 



band to intra-band excitations now occurs at a higher 
interaction strength of I = 3.0. 

Having examined the nature of the eigenstates, we now 
turn our attention to the eigenvalues, and in particular, 
to the neutral excitation energies. For the equal vf sys- 
tem, the ground state energy, first neutral excited state 
energy and the excitation energy for various inte raction 
strengths are calculated and shown in figure 6(a). In the 
non-interacting limit, the excitation energy AE is equal 
to 2tx, as would be expected for a transition between 
bonding and anti-bonding bands. As I is increased, AE 
also increases, until it becomes saturated at 7 as 2.0 at a 
value of 2itvf/L. This corresponds to the discontinuity 
in the charge distribution at 7 w 2 (figure 5(a)) already 
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associated with the change in nature of the excitation 
from inter- to intra-band. 

For the system with Vp < Vp, however, the noninter- 
acting AE is no longer 2t± but is approximately 0.81, 
as indicated by the band diagram (figure ^) . As / is in- 
creased, AE decreases and approaches zero in the range 
0-5 ^$ I S 1-0 j the same region as is dominated by the 
covalent states in the charge distribution. For I > 1.0, 
AE once more increases until it too becomes saturated at 
/ ~ 3.0, at a value of around 1.6. This is not, as might be 
suggested by the fact that all the charge is on the non- 
interacting chain, 2nvp/L, but is in fact slightly more 
than the value of 27Tfj^ ondmg /L. The excitation energy 
also continues growing for / > 3.0, albeit very slowly: 
this, and the fact that AE is not exactly 2irv F ° ndms / L at 
saturation, is because of the combined effects of increase 
in energy due to induced boson states on chain B and 
increase in energy due to chain A interactions overtaking 
the effect of charge removal from chain A to chain B. 



D. Spectral functions 

The spectral functions for the model are now consid- 
ered. Results for the calculation of the spectral func- 
tion for both the equal vf and the Vp < Vp systems 
are presented in figure ^. For the equal Vp system, fig- 
ure |7(a)| , there are three sources of broadening of the 
peaks away from the delta function we would expect for 
a Fermi liquid. The first of these is owing to the finite 
value of the imaginary energy r\ used in the numerical 
calculations, and leads to broadening of the peaks into 
Lorentzians with width ~ r\. The second results from 
the splitting into bonding and antibonding bands due to 
the interchain coupling and is equal to 2tj_. However, the 
remaining broadening, which is interaction-dependent, is 
due to the removal of spectral weight from the Fermi en- 
ergy at w = 0, and is indicative of residual Luttinger 
liquid behaviour in the coupled sy stem . 

For the vp~ < Vp system, figure 7(b) , we again see the 
broadening of individual peaks into Lorentzians. How- 
ever, the peaks initially converge as / increases, analo- 
gous to the similar decrease in neutral excitation energy 
observed in the same region of / < 0.5. In this region 
around / ~ 0.5 however, two "branches" of the spectral 
function emerge, corresponding to separate peaks aris- 
ing from covalent basis states and ionic chain B basis 
states — at lower values of I the covalent peak dominates, 
whereas for higher values of / the ionic B states become 
dominant. For the equal vf spectral function, a similar 
analysis reveals that the single peak is due to ionic B 
states and there is no corresponding second peak aris- 
ing from covalent states. The peak separation at I = 
is approximately 0.81, consistent with that indicated by 
the band structure, and following the emergence of the 
two branches noted above, the peaks separate in the 
I > 1.0 region, again indicating residual Luttinger liq- 
uid behaviour. 
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FIG. 6: Eigenvalues for the ground state and neutral first 
excited state of the two-electron, two-chain model. 



E. Calculation of the Luttinger parameters 

We can use the spectral function results to calculate 
the Luttinger parameters of the system. For a contin- 
uum system, this could be achieved by direct numerical 
measurement of the exponent a |28|, [3(J, but this is 
not possible in small discrete systems such as those under 
present consideration. However, we can define an effec- 
tive Luttinger velocity Vq for the system as follows. For 
a single Luttinger liquid chain, the separation Alu of the 
electron and hole contributions to the spectral function is 
equal to 2v q, whereas for a non-interacting do uble c hain 
system Alu is equal to 2t± . We see from figure 7(a) that 
for the equal vf system, we have a combination of both 
these effects, and hence define our effective vq in terms 
of the electron-hole peak separation 



Alu 



2{vf + tj 



)• 



(24) 
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FIG. 7: Spectral functions p(q = 0, u) for both types of two 
chain system. The hole and electron contributions have been 
renormalized to the same height in order to show the positions 
of the peaks accurately. Creation and annihilation both took 
place on chain A. 



As we have chosen gziq) to be equal to 54(g), we can 
equate equations @ and ([lO]) to give our effective Lut- 
tinger parameter 



Kf 



.20 



VF 
„,cft" 



(25) 



The results of this are shown in figure [fC] and compared 
to the single-chain continuum limit. 



system the situation 



However, for the 
is slightly different, fn the non-interacting limit the 
bonding-antibonding splitting is no longer 2t^, but in- 



stead is of the order of q\vp 
case, much larger than 2ij_ 



■ Vp\. Since this is, in this 
the bonding/antibonding 



splitting may be ignored, resulting in the following ex- 
pression for the effective Luttinger velocity 




B 



(26) 



We may then proceed as before to calculate Vq S and 
hence K p — v s p STEM /vQ S for this system. 



electron contribution 
hole contribution 




FIG. 8: Spectral function p(q — ty/L,uj) for the Vp = 1.0, 
vp — 2.0 two chain system. The hole and electron contribu- 
tions have been renormalized to the same height in order to 
show the positions of the peaks accurately. 



However, we note from figure |8j that there are two sep- 
arate peaks in the electron part of p(q,tu). A decom- 
position of the basis states contributing to these states 
reveals that the lower energy peak is due to chain B-type 
ionic states, while the peak at higher u> is due to the cova- 
lent states, with virtually no contribution at all from the 
chain A ionic states. This is shown in figure ||. At lower 
interaction strengths, the covalent states dominate, but 
as the interaction strength increases, the chain B ionic 
state grow in importance and the ionic states decrease, 
fading completely by / w 2.0. The same phenomenon 
is responsible for the double peak structure in the hole 
contribution to p{q = 0, w), where the effect is somewhat 
less pronounced. 

Either of these peaks can be used to calculate K p . Fig- 



ure 10 shows the results for both the covalent and ionic 
peaks, together with the calculated Luttinger parameters 
for the L = 4 and L = 6 equal chain systems. The values 
for the continuum limit of a single chain are also shown 
for comparison. We can see that for / < 0.5, this system 
is more Luttinger liquid- like than the equal Dp system, 
but for interaction strengths greater than / > 1.0, the 
covalent branch peters out and the system becomes less 
Luttinger liquid-like than the corresponding equal chain 
system. 
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chain B ionic states - — full electron contribution 1/L 



FIG. 9: Decomposition of the electron contribution to the 
spectral function into parts due to ionic chain A states, ionic 
chain B states, and covalent states. For the L = 4, N = 2, 
t± = 0.1, v£ = 1.0, «# = 2.0 system at I = 0.7. 



FIG. 11: Scaling with respect to inverse system length for 
the calculated values of K p for the Vp = Vp = 1.0 system 
at various interaction strengths. The 1/L — data are for a 
single-chain continuum system. 
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FIG. 10: Calculated values of the Luttinger parameter K p for 
several systems. See also reference p2j]. 



IV. CONCLUSIONS 

The basic question we have asked in the current work is 
whether the coupling of an Luttinger liquid wire to a sur- 
face destroys any or all of its one-dimensional properties. 
We have attempted to answer this question purely by 
examination of the eigenstates and correlation functions 
of our, albeit simplified, model system. A more direct 
way of answering this question would be to measure the 
transport properties directly. However this would require 
substantial changes to the model which we believe would 
run the risk of obscuring the basic question: whether 
one-dimensional metallic systems coupled to surfaces in 
general retain their one-dimensional properties, regard- 
less of the influence of such external sources as leads. 



We have presented results for a system of two chains, 
identical except for the presence of electron-electron in- 
teractions on one chain, representing the wire. We have 
shown that the coupling to the surface is weak in na- 
ture even for coupling parameters as high as t± = 0.5, 
as the contribution to the ground state from basis states 
involving bosons remains comparable to that for an iso- 
lated Luttinger-liquid chain. We have seen that there is a 
fundamental change in the nature of the low- lying neutral 
excitations at a certain value of the interaction strength, 
owing mainly to the finite length of the chains. However, 
this shift, from inter-band excitations with Afe = to 
intra-band excitations with Ak = 1~kvf/L also leads to a 
transfer of charge away from the interacting "wire" chain 
to the non-interacting "substrate" chain, resulting in the 
restoration of some Fermi liquid properties. In particular 
the excitation energy in this region depends not on the 
Luttinger velocity vq but on the Fermi velocity vf ■ How- 
ever, the spectral properties, which depend on the nature 
of the charged rather than the neutral excitations, indi- 
cate that Luttinger liquid behaviour survives at values 
of the interaction strength not only below this threshold, 
but also above it. This behaviour manifests itself in the 
removal of spectral weight from the Fermi surface, and 
we have used this feature to calculate effective values of 
the Luttinger parameter K p for our system. 

Our calculated values of K p indicate definitive Lut- 
tinger liquid behaviour for the system as a whole, and 
moreover, tend toward the values for an isolated contin- 
uum chain as the system length L increases. Figure |ll| 
shows the scaling of K p as a function of inverse system 
length for the equal- «f system at various values of the 
interaction strength /. The values for K p for an isolated 
continuum Luttinger liquid (L — 0) are included and a 
straight line fit is shown through all three data points. 
One would not expect a perfect straight line fit from these 
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data — the continuum single chain system remains, after 
all, substantially different in nature from a continuum 
two-chain system. However, the fact that the data are 
so linear with respect to inverse system length further 
reinforces our general conclusions. 

In summary, therefore, we have investigated the prob- 
lem of electron-electron interactions in an atomic scale 
wire on a surface. We have presented results that indi- 
cate that the resulting system does not behave as a pure 
Luttinger liquid but nonetheless retains many character- 
istics consistent with Luttinger liquid behaviour. It is 
clear from our results, however, that there is an impor- 
tant caveat to be borne in mind for hybrid systems such 
as this — that the results obtained depend very much on 
the properties of the system one chooses to measure. 
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